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Summary 

A new hybrid surface-integral-finite-element numerical scheme has been developed to model a 
three-dimensional crack propagating through a thin, multi-layered coating. The finite element method 
was used to model the physical state of the coating (far field), and the surface integral method was used 
to model the fatigue crack growth. The two formulations are coupled through the need to satisfy 
boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak 
that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume 
can be set up independently. Thus when modeling crack growth, the finite element mesh can remain 
fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented 
to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of 
surface cracks propagating in engine components. In this work, the authors formulate the hybrid 
surface-integral-finite-element method and discuss the mechanical issues of implementing a structural 
health monitoring system in an aircraft engine environment. 

Introduction 

Most field failure issues of components can be traced back to design, manufacturing, repair and 
inspection processes or unidentified usage conditions. In severe cases, these conditions can cause 
fatigue cracks to initiate and propagate, sometimes undetected, through the part causing catastrophic 
failure. Fortunately, most of these cracks are detected through routine inspection and overhaul. 
However, there are instances where abnormalities in the manufacturing process can cause undetectable 
imperfections that may produce fatigue cracks in the fielded part. Therefore, use of structural health 
monitoring systems capable of detecting impending failures would allow timely maintenance actions to 
increase overall safety. The economic impact of such a system is tremendous, based on four principal 
effects: (i) reducing the inventory, (ii) reducing maintenance cost, (iii) extending life of structures by 
condition based maintenance, and (iv) market advantage. 

In this paper, the authors develop an effective numerical procedure to solve spatial problems of 
crack propagation in heterogeneous media and propose a structural health monitoring system. Work 
focuses on the development of a hybrid numerical method to analyze a crack propagating through a 
thin, multi-layer coated part. The primary objective of the hybrid method is to combine the advantages 
of finite elements and boundary elements to predict three-dimensional crack propagation. The finite 
element method (FEM) was used to model the state of the part (far field), and the surface integral 
method (SIM) was used to model the fatigue crack growth. A key advantage of this approach is that 
only the crack surfaces have to be re -meshed during crack propagation, with the FEM remaining the 
same. 

The authors further propose a structural health monitoring system that operates by detecting a 
change in the electrical impedance of a conductive mesh due to mechanical changes caused by cracks. 
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The model crack sensor consists of three layers deposited sequentially onto the surface of the 
component, as shown in Figure 1. The first insulating layer is deposited directly onto the component to 
electrically isolate a conductive mesh from the part. The second layer consists of a metallic wire grid 
deposited over the insulating layer. Finally, a protective topcoat is deposited over the metallic grid to 
protect it from oxidative or chemical attack, as well as to protect the grid from electrical interference 
due to conductive particle contamination. Cracks propagating at the component surface will fracture 
the grid resulting in impedance changes of the grid that will be detected by a sensor electronics 
package. The spatial separation between wires determines the crack size resolution the mesh can 
detect. 



Figure 1 : 3D problem formulation - geometry and 
loading conditions. 


Formulation of the Numerical Method 

The surface integral method is an indirect boundary element method based on the superposition of 
force multipoles to represent two-[l, 2] and three-dimensional [3, 4, 5] fractures. The tractions induced 
by a fracture can thus be represented by the integral equation 

tk = \\/bnn n m 8 n dA ( 1 ) 

where t k is the traction vector, A is the surface area of the nonplanar fracture, r kmn is the stress influence 
function, n m is the normal to the fracture surface, and 5„ represents variations of crack displacement 
where is crack opening, and S 2 , S 3 are orthogonal components of crack shear. The stress influence 
function used in this formulation was obtained by differentiating Rongved's fundamental solution [6] 
for point forces near a planar bimaterial interface. 

Traction boundary conditions are enforced at collocation points located at the element centroids of 
a piecewise planar representation of the fracture. In the vicinity of a collocation point, one can subtract 


off the integral equivalent of a rigid body displacement to reduce the order of the singularity that can be 
integrated in the Cauchy principal-value sense. Discretization of the fracture surface into crack 
elements leads to the following equivalent form of the traction equation 
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where Nele is the number of elements, the superscript i indicates crack element correspondence, q 
identifies the singular crack element containing the collocation point, h is the local interpolation 
function defined such that the near-tip displacements vary as the square root of the perpendicular 
distance from the crack front, A t is the surface area of the z'th crack element, N is the outward unit 

normal to the contour enclosing A q , and it follows from Green's Lemma that ^ kmn n q m = V ■ V0 A .„ where 

r knm and § kn are continuously differentiable along the closed contour. Using shorthand notation, 
equation (2) can be rewritten as 

W=H<5} (3) 


where C kn represents the coefficient matrix. A full derivation of the governing equations can be found 
in Forth and Keat (1996) [4], 

The combination of the surface integral method with the finite element method has been 
accomplished for two- [2] and three-dimensional [7] fracture analyses using superposition. The two 
constitutive models are coupled through the need to satisfy boundary conditions on the external 
boundary and crack surface. The formulation of the finite element equations is well documented, and 
can be presented in shorthand notation as 

rfy ra }={i?}-{tf c } (4) 

where [ K\ is the finite element stiffness matrix, {U FE } is the vector of finite element nodal 
displacements, {R} is the external nodal load vector, and {R c } is introduced to the equation to satisfy 
the external traction boundary conditions that are not enforced in the surface integral model. {R c } can 
be defined with respect to the surface integral crack opening displacement vector such that 

(5) 

where [G] is derived such that the nodal forces are statically equivalent to the surface integral tractions 
acting on the external surface patches. 


The traction equation for the surface integral method, as defined in equation (3), then can be 
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modified to accommodate the finite element correction {r } at the crack. 

{ t c } can then be expressed in terms of the finite element nodal displacements such that 

{t c }=[5]{(7 F£ } 
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The complete system of equations can be written by combining equations (4), (5), (6) and (7) to obtain 
the following partitioned matrix equation 
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To impose displacement boundary conditions, the equation must be reformulated in tenns of the total 
displacements { U TOT } as 

{t/ ror }={E/™}+ {[/*'} 


where 
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in which [X] is computed with the displacement influence functions. Therefore, the final form of the 
hybrid governing integral equation can be written as 
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Crack Propagation Model 

Three-dimensional fracture models, such as the surface integral method, can make use of two- 
dimensional growth laws by applying them locally along the crack front. To propagate the crack, both 
the direction and extension must be computed at each crack tip element. The maximum circumferential 
stress theory [8] for direction, and the Fonnan fatigue crack growth equation [9], were applied here due 
to their wide industrial acceptance and ease of implementation. Maximum circumferential stress theory 
determines the direction of crack propagation, 0, by maximizing the circumferential stress at the crack 
tip such that: 


Kj sinO + K H (3cos# - 1) = 0 
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where 0 is the angle of growth direction expressed in the local nonnal-tangential coordinate system of 
the tip element and K] and K n are the mode I and mode II stress intensity factors, respectively. The 
Fonnan equation was developed to curve fit a large material database that is available through NASA 
Johnson Space Center and is defined as: 
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where da/dN is the rate of crack growth based on crack length a and cycle count N; R is the stress ratio 
(5 mi max 5 AAA is the equivalent stress intensity factor range which for mixed-mode is defined to be 


AK eq = ^(AKj ) 2 + {AK U ) 2 + (A Kjjj f in accordance with the mixed-mode definition of strain energy 
release rate [10]; C, n, p and q are empirically derived material constants; K c is the fracture toughness; 
AK th is the threshold stress intensity factor range and /is the crack opening function. 




Figure 2. Hybrid Finite element mesh (FEM) and Surface integral mesh (SIM) to predict 
wired coating failure. Wires break as the crack in the part propagates. 



Figure 3. Hybrid Finite element mesh (FEM) and Surface integral mesh (SIM) to predict 
wired coating failure. Fracture of the coating system due to thickness and wire spacing. 


Structural Health Monitoring System Analysis 

The numerical model was comprised of 240 3D brick, 8-noded finite elements and 200 surface 
integral elements and the crack was propagated using material properties and loading conditions 
representative of an engine operating environment. The initial crack radius was defined as a 0 , the base 
material (E b rq a! are Young’s modulus, Poisson’s ratio and coefficient of thermal expansion, 
respectively) and the wire mesh were assumed to be of the same material. In addition, the insulating 
and protective layers were assumed to have the same material properties (E 2 , V 2 , 0t 2 ) The initial crack 
was located equi-distant between wires such that the wire spacing (s) would be minimized to detect a 
crack of size a 0 as it propagated. All spatial variables are non-dimensionalized to the initial crack 
length, do, and material variables are non-dimensionalized to the component material (E b \q oq). 
Initially, the crack has a penny shape and is completely buried in the base material. 


In this paper the authors numerically studied two crack propagation inodes, specifically through- 
layer cracking (Fig. 2) and coating failure (Fig. 3). The insulating layer and protective top-coat are 
brittle, e.g. a lower fracture toughness, in comparison to the component material and will therefore 
fracture at a lower stress. The key factors in preventing fracture of the coating are the position and 
spacing of the wire in the mesh. Further details of this structural health monitoring system can be 
found in Forth and Staroselsky (2001) [11]. If the wire is placed too close to the surface, the insulating 
layer will simply fall off as the fatigue crack propagates through the part, while the wire remains intact. 
If the wire is placed too close to the part, then the fatigue crack will break through the insulating layer 
between the part and the wire causing delamination to occur at the interface. If the wires are spaced too 
far apart, a crack of size a 0 will not be detected. Therefore, the wire in the sensor must be embedded 
within the coating a minimum of 25% of the coating thickness from either surface to provide adequate 
structural integrity and electrical insulation between the component and wire. The wires must also be 
spaced a maximum of eighty percent of the predetermined detectable crack size, a 0 . 

Conclusions 

A novel hybrid numerical method to study three-dimensional crack propagation in a layered 
structure has been developed. This method combined with a crack propagation criterion made possible 
the study of incremental growth of three-dimensional non-planar cracks as well as to hypothesize the 
requirements for a structural health monitoring system. The fracture mechanics modeling has shown 
that it is feasible to design a patterned coating on an engine component which will be adherent at the 
stresses and temperatures normally present in engine operation, and will fracture in a controlled manner 
when there is an underlying crack in the part. The model has predicted an upper boundary for the total 
coating thickness, of which at least 25% of the coating thickness must be the protective layer covering 
the conductive grid, and a wire spacing of 80% of the predetermined detectable flaw size for the system 
to be structurally reliable. 
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